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DESCRIPTION OF A COMPUTERIZED METHOD FOR PREDICTING 
THERMAL FATIGUE LIFE OF METALS 
by D A, Spera and E C. Cox 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 4^135 

ABSTRACT 

A computer program called "TFLIFE" is described which can be used 
to predict the thermal fatigue life of metals and structural components 
from conventional metal properties. This program is used as a 
subroutine with a main program supplied b> the user. The main program 
calculates input cycles of temperature and total strain for TFLIFE which 
then calculates a stress cycle, creep and plastic strain damage, and 
cyclic life. A unique feature cf TFLIFE is that it calculates lives 
according to several different failure criteria for the same input data 
These criteria are surface crack init lar ion, interior crack initiation, 
and complete fracture of both unnotched and notched fatigue specimens. 
Sample output tables ate shown, together with results for two typical 
problems: (1) thermal-mechanical latigue oi bar specimens of the 

tantalum alloy T-lIl, and (2) thermal-stress fatigue of wedge specimens 
of the nickel alloy B-I900. Thermal latigue lives calculated using 
TFLIFE have been verified by comparison with a variety of laboratory 
test data on different types ot alloys The computer program is now 
ready for mote extensive evaluation on structural components as well as 
additional laboratory specimens. 


INTRODUCTION 

This paper describes a fatigue life calculation system utlll^lng a 
computer program called "TFLIFE". TFLIFE calculates the low-cycle 
fatigue life of metals aubjerted to simultaneous thermal and strain 
cycling. In conjunction with an elastic structural analysis program, 
TFLIFE can be used to predict the fatigue life of complete components in 
high-temperature equipment such as engines and heat exchangers This 
paper contains a general description ct the program, typical output 
tables, and results for two sample laboratory test problems. 

The original version of TFLIFE was written more than eight years 
ago to aid in the development of a .teep damage theory for thermal 
fatigue (refs 1 and 2) Since 'hat time it has grown in scope and 
function, serving as a useful tool tor research on high-temperature 
metals (ref 3). It now in. lodes recent developments in creep-fatigue 
theory and input/output formats which make it suitable for general use 
To date, TFLIFE results have been verified by thermal fatigue tests on 
laboratory specimens such as oars, tapered disks, wedges, and simulated 
turbine blades (rets 1 - 4 ) In these tests thermal and mechanical loads 
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were a'^plled by fatigue machines, fluidized beds, and burner rigs to 
specimens of nickel. Iron, cobalt, and tantalum alloys. TFLIFE Is now 
ready for more extensive evaluation by the technical community on 
complete components as well as laboratory specimens. 

A unique feature of TFLIFE Is that It calculates thermal fatigue 
lives according to each of three failure criteria, based on one set of 
input data. These criteria are (1) surface crack Initiation, (2) 
interior crack initiation (applicable to coated metals) and (3) complete 
fracture of thermal-mechanical fatigue specimens, both unnotched and 
notched. Including these three criteria in a single computer program 
produces a very general and complete analysis of thermal fatigue 
resistance . 

A complete set of the equations used in TFLIFE is given in Ref. 5, 
together with derivations. Computer time for a typical problem is 8 to 
17 seconds on an IBM 360/67 time-sharing system. Required storage 
capacity is approximately 12000 words. 

GENERAL DESCRIPTION OF TFLIFE 

Figure 1 is a schematic diagram which describes the general 
sequence of operations performed by the life calculation system. TFLIFE 
is a subroutine which is called by a user-supplied main program. This 
main program, discussed below, contains all calculations of temperature 
and total strain specific to the specimen or component being analyzed. 
Initially, the main program calls TFLIFE which reads its own input data 
(Block I) and divides both the heating and cooling periods into 
approximately 100 Increments. The input data for TFLIFE are of two 
types: (1) conventional mechanical properties which describe the 

monotonic tensile and creep-rupture behavior of the metal, and (2) 
optional empirical constants based on cyclic data. The conventional 
metal properties must be given for a range of temperatures which 
Includes the therma. cycle being analyzed. The optional empirical 
constants pertain to cyclic strain hardening, strain concentration, and 
time- Independent fatigue behavior. Nominal values for these constants 
are contained in TFLIFE for cases in which cyclic data are not 
available. Input data can be in any of the following three systems of 
units for stress and temperature: (1) pounds per square inch and 

Fahrenheit degrees, (2) grams per square millimeter and Celsius degrees, 
or (3) Newtons per square centimeter and Kelvin degrees. SI (Systems 
International) units can be specified for the bulk of the output, 
irrespective of the input unit system. TFLIFE then returns to the main 
program. 


Temperature, Strain, and Stress Cycles 

At Block II in the main program, cycles of temperature and 
mechanical atrain (total strain less free thermal expansion) must be 
calculated for a potential failure location in the specimen or 
structure. The complexity of these computations can vary greatly, 
depending on the structure being analyzed. For a strain-cycled 
thermal-mechanical fatigue specimen. Block II may only be required to 
generate saw-tooth or sinusoidal cycles. For a thermal-stress fatigue 


specimen it may be necessary to calculate transient temperature 
distributions and the strains which result from them. In this case a 
linear elastic analysis is usually sufficient, according to the 
hypothesis of total strain invariance proposed in Ref 6. This 
hypothesis states that under conditions of thermal stress, the 
displacements in a body are substantially independent of local inelastic 
strains. Thus, displacements obtained by a linear elastic analysis can 
be used to calculate mechanical strains which can later be divided into 
elastic, plastic, and creep components. For a complex component it may 
be necessary for the main program to call a finite-element analysis 
subroutine. In all cases, a local temperature and a local uniaxial 
mechanical strain are calculated for each time increment. The system 
then calls TFLIFE a second time. 

In Block III an incremental stress-strain analysis is performed 
which divides the mechanical strain cycle into elastic, plastic, and 
creep components. This analysis is repeated until "shakedown" occurs, 
when the stress-strain hysteresis loop becomes stable and repetitive. 

In this way a stress, a plastic strain increment, and a creep strain 
increment are calculated for each time increment. Approximately 
two-thirds of TFLIFE's computing time is used in these calculations. 

Upon completion of Block 111 the mechanical behavior of the material at 
the potential failure location has been thoroughly described in terms of 
temperature, stress, and strain Fatigue damage can now be calculated. 

Fatigue Damage and Life Predictions 

In Block IV TFLIFE calculates the fraction of the metal's life 
consumed in one cycle This calculation uses a fatigue damage model 
which we call the generalized damage fraction theory (ref 5). The model 
is so named because it integrates elements from the work of many 
investigators into a unified, practical method of life analysis. This 
theory uses Palmgren's linear '■emulative damage hypothesis (ref 7) as 
modified by Taira to include time as a variable (ref 8). Fatigue damage 
is assumed to to be ot two types: (1) elastic-plastic strain damage 

which is independent of time, and (2) creep strain damage which is 
t ime-dependent 

Elastic-plastic strain damage is calculated using an empirical 
fatigue equation of the exponential type, derived from equations 
proposed by Basquin (ret 9 ), Manson (refs 10,11), and Coffin (ref 12). 
Equations for creep strain damage combine contributions by Robinson (ref 
13), Hoffman (ref 14), Taira (ref 8), Swindeman (ref 15), and Spera 
(refs 1 to 5) . 

In Block V life predictions are now made for each of the following 
three tail ire criteria: (1) surface crack initiation, (2) interior 

crack initiation, and (3) complete fracture. The first criterion 
applies mainly to thermal-stress fatigue of an uncoated material. The 
second criterion can be used to predict the thermal-stress fatigue life 
of a coated material in which the fatigue life of the coating exceeds 
that ot the substrate material. The third criterion is limited to 
thermal-mechanical tatigue testing in which strains and temperatures are 
nominally uniiorm in the test section of the specimen. The program 
provides tor three strain concentration factors so that the user can 
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calculate the effect of notch severity on fatigue life. One of these 
factors is always unity, representing the unnotched condition. The 
other two factors may be specified arbitrarily by the user or nominal 
values of 2.0 and 3.0 will be supplied by the program For notched 
specimens the input temperature and mechanical strain cycles are 
considered to be nominal rather than local. Thus five cyclic lives are 
calculated in Block V for one pair of temperature and mechanical strain 
cycles as calculated in Block 11. This is a unique feature of the 
TFLIFE computer program. 


Output Tables 

In Block VI the results of the life analysis are printed in the 
form of six tables. Samples of some of these tables are shown in Fig. 2 
for a thermal-mechanical fatigue problem which will be discussed later. 
Table 1 (not shown) documents the applied cycle in terms of the 
mechanical strain range plus temperature, strain, and time limits. 

Table 2 contains the actual life predictions. These are presented for 
each of the three failure criteria discussed above. This table also 
contains damage fractions which show the theoretical proportions of the 
damage associated with creep strain and elastic-plastic strain for each 
predicted life. The sum of the two damage fractions is always unity. 
Also, interior crack initiation and complete fracture are assumed to 
occur simultaneously in a thermal-mechanical fatigue specimen because of 
the nominally uniform strain and temperature conditions in the test 
section. 

Table 3 presents temperature, strain, and stress as functions of 
time during heating and cooling. If the user wishes, this table can be 
doubled or tripled in length, resulting in smaller time increments and 
more detailed data on the cycle. Table 4 lontains the monotonic tensile 
data for the material as a function of temperature This table is 
constructed from the actual input data using linear interpolation The 
number, order, and spacing of temperatures for which property data are 
given are arbitrary. 

The cyclic hardening ratio in the footnote to Table i* relates 
cyclic and monotonic stress-strain data. This ratio is an optional 
scaling factor for proportionately increasing or decreasing yield and 
ultimate tensile stress values at all temperatures. It is specified by 
the user according to one of the following three criteria: (1) If no 

cyclic stress-strain data are available the hardening ratio is unity, 
which represents neutral hardening. Also, unity is the nominal value 
used by the program in the absence of a specified value (2) If partial 
cyclic stress-strain data are available (for example, at only one 
temperature) the hardening ratio is adjusted until flrs^-ctder 
differences between calculated and observed stresses are eliminated. 

(3) If complete cyclic stress-strain data ate available they are used to 
calculate equivalent monotonic data for a cyclic hardening ratio of 
unity. 

Table 5 is a summary of the empirical equations which are used in 
TFLIFE to calculate creep and rupture times (ref 1). The coefficients 
are calculated from test data by means of regression analysis. Most of 
the commonly used time-temperature parameters are included in these 
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general equations as special cases, either exactly or to a close 
approximation. In this way the equations in Table 5 represent an 
attempt to select the optimum time-temperature parameter for the 
available data. It only limited creep rate data are available, 
simplifying assumptions can be made by the user so that some of the 
curve-fit constants in the equation for time to one-percent creep are 
the same as those in the equation for time to rupture. 

A final table. Table 6 (not shown), contains data which document 
the Iterations required to acheive convergence. Each iteration is a new 
calculation of the stress and inelastic strain cycles. The stress at 
the end of the cooling cycle is carried forward to start the incremental 
stress analysis for the next iteration. A residual Inelastic strain is 
recorded which is the sum of the plastic and creep strains at the end of 
each iteration. (Both these quantities are set to zero at the start of 
each Iteration.) A residual inelastic strain of zero represents a 
completely stable, repetitive stress-strain hysteresis loop. When the 
residual strain becomes less than 10 micro-units, a curve-fit procedure 
IS used to accelerate convergence, and iteration Is terminated after two 
additional cycles. 


Options 

After output tables have been printed for one set of input data, 
there is a minimum of three options which the user can choose to Include 
In the main program: (1) terminating the program, (2) starting the 
solution to a new problem by reading a new set of data, or (3) varying 
parameters In the input data which are already stored in the computer 
and starting the solution to a modified problem. The third option 
enables the user to perform a parametric study from a single set of 
Input data. As mentioned previously, the operations In the main program 
are determined by the needs of the user. The TFLIFE system shown 
schematically in Fig. 1 Includes sufficient flexibility to satisfy a 
wide variety of such needs. 


SAMPLE RESULTS 

The application of TFLIFE to different types of thermal fatigue 
conditions will now be illustrated by showing some typical results for 
the following two problems: (1) thermal-mechanical fatigue of bar 

specimens of the tantalum-base alloy T-111 tested in vacuum (ref 16), 
and (2) thermal-stress fatigue of wedge specimens of the nickel-base 
alloy B-1900 tested in fluidized beds (refs 17,18). Figure 3 shows the 
geometry of these two types of specimens. The slot notch shown in Fig. 
3(b) has a theoretical elastic stress (or strain) concentration factor 
of 3.0 (ref 19). Predicted lives are given for unnotched and notched 
bars of T-111 as a function of imposed strain range combined with either 
isothermal or in-phase thermal cycling. The B-1900 vadge specimen 
calculations Include a study of the effects of changes in cycle time and 
edge radius for both the uncoated and alumlnlde-coated conditions of the 
alloy. 
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Thermal-Mechanical Fatigue of T-lll Bars 

i 

' T-lll was selected as a sample allo> because of the variety of 

thermal-mechanical fatigue data reported cn it in Ref 16. Also, its 
I stress-strain behavior Is as complex as any which the user is likely to 

encounter. Therefore, T-lll provides an excellent test of the 
capabilities of TFLIFE with respect to thermal-mechanical fatigue. The 
output tables presented in Fig. 2 are taken from one of the T-lll test 
j cases. 

Figure A shows the simple type of mechanical strain and temperature 
cycles programmed for these tests. When diametral strain is the control 
variable, the normal input variable of axial strain is usually unknown 
prior to the stress analysis (Fig. I, Block III). TFLIFE resolves this 
difficulty by using an equivalent axial strain equal to twice the 
diametral strain plus an equivalent elastic modulus equal to Young's 
modulus divided by twice Poisson's ratio (r^'f 5) 

Figure 5 illustrates the stress-strain behavior calculated for both 
the Isothermal and in-phase thermal cycles with equal strain ranges. 

The significantly different shapes of these two calculated hysteresis 
! loops are confirmed by test data (ref 16). Cyclic hardening ratios were 

selected according to the criteria described previously for cases with 
I partial cyclic stress-strain data (roe Output Tables). Two different 

cyclic hardening ratios were required for these two cycles because the 
stresses in the in-phase thermal cycle were found to be significantly 
larger than those in the isothermal cycle for equal strain ranges. A 
' hardening ratio of 1.43 was needed to correlate calculated and observed 

I stresses at point A in the in-phase loop. A much smaller factor, 1.08, 

matched the calculated stress amplitude in the isothermal loop (point B) 
I to the measured value. These hardening ratios reflect the complex 

cyclic hardening which has been observed in T-lll. For behavior of this 
j type an accurate estimate of the hardening ratio is essential to 

accuracy in the prediction of creep damage and life. 

^ Figure 6 shows the results of the life calculations for isothermal 

I cycling as a function of axial mechanical strain range. The failure 

I criterion for these tests was complete fracture cf the bar specimen. 

Predicted behavior with a strain concentration factor of unity is in 
good agreement with the unnotched specimen data, as shown by the open 
symbols. This correlation between theory and experiment gives a measure 
a confidence in the extrapolated behavior for lives greater than 1000 
cycles. For the notched specimens (closed symbols) predicted behavior 
is somewhat conservative when the strain concentration factor is equated 
to its theoretical elastic value of 3.00. The available data were best 
fit using an empirical strain ccncentraticn factor of 2 15. Life 
predictions based on this factor are given by the dashed line in Fig. 6. 
In the absence of notched specimen data, the theoretical elastic strain 
concentration factor would be used as input in TFLIFE. 

Tn Fig. 7 predicted and observed lives are compared for tests with 
In-phase thermal cycling. Unlike the isothermal tests, the in-phase 
tests of T-lll showed little notch effect. The TFLIFE predictions are 
in agreement with the data, whether the strain concentration factor is 
assumed to be 1.00, 2.15, ct 3.00. This insensitivity to strain 
I concentration is related directly to an assumption In the creep damage 

I model in TFLIFE. The assumption is that interior creep damage is not 
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affected by a surface strain concentration but depends only on the 
nominal stress. Thus, the strain concentration affects only the 
elastic-plastic strain component of damage. During In-phase thermal 
cycling at relatively large strain ranges (lives less than 1000 cycles) 
most of the damage results from creep (Fig. 2, Table 2). Therefore, the 
model Indicates that strain concentration has only a minor effect on the 
predicted lives at these strain ranges for In-phase cycling. However, 
as the strain range decreases below about 0.004, elastic-plastic strain 
damage becomes dominant over creep damage and a significant notch effect 
Is predicted for this alloy. Figure 7 is an example of an extrapolation 
of data by TFLIFE that indicates a notch effect at long life which is 
not present In the short-life data. 

Thermal-Stress Fatigue of B-1900 Wedges 

Rapid heating and cooling of the wedge specimens shown in Fig. 3 
produces edge cracks similar to those often found in turbine blades and 
vanes. The sequenc-j of calculations required for predicting the 
Initiation period for such cracks is illustrated by the results shown In 
Figs. 8 to 11. In contrast to the simple cycles shown in Fig. 4, the 
transient temperature and strain inputs for thermal-stress fatigue 
problems are quite complex. Consequently, Block II in the main program 
(Fig. 1) Is much more complex for thermal-stress fatigue than for 
thermal-mechanical fat...gue problems. For Instance, in this B-1900 
problem, thermocouple data taken at several points In the interior of 
the wedge specimen are curve-fit in Block II for each of the 
approximately 100 time Increments during heating, and the same number of 
time increments during cooling, as described in Ref 16. These 
curve-fits are then used to calculate the temperature and mechanical 
strain cycles at the specimen edges. 

Figure 8 shows typical edge temperature transients during heating 
and cooling. Note that the length of the heating and cooling periods Is 
not sufficient for the temperatures to reach steady-state levels equal 
to the bed temperatures. In Fig. 9 transient mechanical strain is 
plotted versus edge temperature rather than versus time. This is often 
a convenient way to display transient data and to insure that calculated 
strain and temperature cycles are compatible. Changes during one of the 
transfer periods are also noted in Figs. 6 and 9, when the specimen Is 
moved from the cooling to the heating bed. 

A typical stress-strain hysteresis loop for thermal-stress fatigue 
is presented In Fig. 10. A cyclic hardening ratio of unity (i.e. 
neutral hardening) was used to calculate this loop, in the absence of 
pertinent cyclic stress-strain data. Comparison of Figs. 5 and 10 shows 
the wide variety of stress-strain behavior possible during thermal 
fatigue. The hysteresis loop in Fig. 10 has a "figure eight" shape 
which Is produced by changes in elastic modulus with temperature 
combined with the presence of only small amounts of inelastic strain. 

The shakedown process during iteration is shown to produce an Increase 
in this cycle's mean stress. 

Figure 11 summarizes the dependence of the number of cycles 
required to initiate edge cracks on the parameters of cycle time and 
edge radius, for uncoated and aluminide-coated conditions. Comparison 
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of predicted curves with test data indicate that TFLIFE can correctly 
model thermal-stress fatigue life on the basis of theimocouple date and 
monotonic material properties. Furthermore, TFLIFE can be used to 
predict behavior outside the range of the test parameters. 

CONCLUDING REMARKS 

The TFLIFE computer program has been found to be a useful tool for 
analyzing thermal fatigue data in terms of conventional metal properties 
and for predicting the thermal fatigue life of metals. It contains 
models for the major damaging phenomena which occur during a general 
thermal-fatigue cycle, at least to a first-order approximation. Thermal 
fatigue lives calculated using TFLIFE have been verified by comparison 
with a variety of laboratory test data on different types of alloys. 

This computer program is ready for more extensive evaluation by the 
technical community, not only on additional laboratory specimens but 
also on structural components. 
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